Achieving a Strongly Temperature-Dependent Casimir Effect 
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We propose a method of achieving large temperature T sensitivity in the Casimir force that in- 
volves measuring the stable separation between dielectric objects immersed in a fluid. We study the 
Casimir force between slabs and spheres using realistic material models, and find large > 2nm/K 
variations in their stable separations (hundreds of nanometers) near room temperature. In addi- 
tion, we analyze the effects of Brownian motion on suspended objects, and show that the average 
separation is also sensitive to changes in T. Finally, this approach also leads to rich qualitative 
phenomena, such as irreversible transitions, from suspension to stiction, as T is varied. 

PACS numbers: 



Casimir forces between macroscopic objects arise from 
thermodynamic electromangetic fluctuations, which per- 
sist even in the limit of zero temperature due to quantum- 
mechanical effects (the Bose-Einstein distribution of the 
photon fluctuations) pQ. In most vacuum-separated ge- 
ometries, such as parallel metal plates, the force is at- 
tractive and decaying as a function of plate-plate sep- 
aration pQ, becoming readily observable at micron and 
submicron separations. For a nonzero temperature T, 
the force is predicted to change as a consequence of the 
changing photon thermal distribution, but this change 
is typically negligible near room temperature and sub- 
micron separations [2j [3] and is only a few percent for 
~ 100 K changes in T at 1-2 jiva separations where 
Casimir forces are barely observable [2-4 . Therefore, 
despite theoretical interest in these T effects [2j [5], it 
has proven difficult for experiments to unambiguously 
observe T corrections to the Casimir force [4]. Other at- 
tempts to measure T Casimir corrections have focused on 
nonequilibrium situations that differ conceptually from 
forces due purely to equilibrium fluctuations [6 . A clear 
experimental verification of a T Casimir correction would 
be important in order to further validate the foundation 
of Lifshitz theory for Casimir effects [2j [3] . 

In this letter, we propose a method for obtaining 
strongly temperature-dependent Casimir effects by ex- 
ploiting geometries involving fluid-separated dielectric 
objects (with separations in the hundreds of nanometers). 
In fluid-separated geometries, the Casimir force can be 
repulsive [HE], an d can even lead to stable suspensions of 
objects due to force-sign transitions from material disper- 
sion [8, 9 or gravity [9j[l0]. We show that, by a proper 
choice of materials/geometries, this stable separation d 
can depend dramatically on T (2 nm/K is easily obtain- 
able), and there can even be transitions where d jumps 
discontinously at some T. Essentially, a stable separa- 
tion arises from a delicate cancellation of attractive and 
repulsive contributions to the force from fluctuations at 
different frequencies, and this cancellation is easily al- 



tered or upset by the T corrections. This appears to be 
the first prediction of a strong T-dependent Casimir phe- 
nomenon at submicron separations where Casimir effects 
are most easily observed. We present the phenomenon 
in simple parallel-plate geometries, but we believe that 
the basic idea should extend to many other geometries 
and materials combinations that have yet to be explored. 
Finally, we also point out that the same systems that are 
strongly T-dependent can also be very sensitive to the 
precise details of the material dispersion at low frequen- 
cies, a property that we plan to exploit in the future. 

The Casimir force between two bodies is a combination 
of fluctuations at all frequencies cj, and at T = can be 
expressed as an integral F(0) = f(£,)d{; over imagi- 
nary frequencies uj = i£ pQ. The contributions /(£) from 
each imaginary frequency are a complicated function of 
the geometry and materials, but they can be computed 
in a variety of ways, such as mean stress tensors and 
the fluctuation-dissipation theorem [7J (valid in fluids for 
subtle reasons [TTJIH]) or via the Casimir energy [12]. At 
a finite T, this integral is replaced by a sum over "Mat- 
subara frequencies" 2i\nkT jh = n£r for integers n: 
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Physically, this arises as a consequence of the 
coth(fkj / '2kT) Bose-Einstein distribution of fluctuations 
at real frequencies — when one performs a contour inte- 
gration in the upper- half complex-cj plane, the residues 
of the coth poles at hu/2kT = nin lead to the sum- 
mation [4]. Mathematically, Eq. corresponds exactly 
(including the x /2 factor for the zero- frequency contri- 
bution) to a trapezoidal-rule approximation to the F(0) 
integral, which allows one to use the well-known conver- 
gence properties of the trapezoidal rule [13] to under- 
stand the magnitude of the T correction. In particular, 
the difference between the trapezoidal rule and the exact 
integral scales as 0(T 2 ) for smooth /(£) with nonzero 
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derivative /'(O) [2] (typical for Casimir forces between 
metals [3]). More specifically, /(£) is exponentially de- 
caying with a decay-length 2irc/a for some characteristic 
lengthscale a (e.g. the separation or size of the partic- 
ipating objects), while the discrete sum of Eq. cor- 
responds to a lengthscale given by the Matsubara wave- 
length At = 2ttc/^t = hc/kT, in which case one would 
expect the T correction to scale as 0(A T /a 2 ). Unfortu- 
nately, at T = 300 K, At ~ 7.6 /im, which is why the T 
corrections are typically so small unless a > 1 /im [2 . We 
cannot change the smoothness of f(uS) since it arises from 
the analyticity of the classical electromagnetic Green's 
function in the upper-half complex-cj plane pQ, so the 
only way to obtain a larger T correction is to introduce 
a longer lengthscale A into the problem that dominates 
over other lengthscales such as the separation a. One 
way of achieving this is to make the /(£) integrand os- 
cillatory with an oscillation period A£ ~ 2ttc/A that is 
much shorter than the decay length ~ 2irc/a. Intuitively, 
discretizing an oscillatory integral induces much larger 
discretization effects than for a non-oscillatory integral, 
and this intuition can be formallized by a Fourier anal- 
ysis of the convergence rate of the trapezoidal rule [15] . 
The question then becomes: how does one obtain an os- 
cillatory Casimir integral? 

One way to obtain oscillatory frequency contributions 
to the Casimir force is to employ a system where there are 
combinations of attractive and repulsive contributions. 
In particular, it is well known that the sign of /(£) be- 
tween two dielectric objects emebedded in a fluid depends 
on the ordering of their dielectric functions at £ [7J [14] : 



Matsubara temperature T (K) 



sgn(/(0) = 



-1, £i(i£) < Efi uid {i£) < s 2 {i£) 
1, otherwise, 
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where +/— denotes an attractive/repulsive force. Since 
the Casimir force depends on the dielectric response of 
the participating objects over a wide range of £, from 
£ = all the way to £ ~ 2itc/a (where a is a characteristic 
lengthscale), the sign and magnitude of the total force at 
any given separation can be changed by a proper choice of 
material dispersion, leading to the possibility of obtain- 
ing Casimir equilibria between objects at multiple sepa- 
rations. This idea was recently exploited to demonstrate 
the possibility of obtaining stable nontouching configu- 
rations of dielectric objects amenable to experiments [9]. 
In this paper, for the purpose of achieving a strong In- 
dependence at short (submicron) separations, we search 
for materials or geometries with dielectric crossings ocur- 
ring at sufficiently small £ = 2ttc/A ~ £t 5 close to the 
room-temperature Matsubara- frequency scale £t- 

To begin with, we compute the Casimir force between 
semi-infinite slabs, computed via a generalization of the 
Lifshitz formula [19] that can handle multi-layer dielec- 
tric objects, with relative permittivities e plotted in Fig.[l] 
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FIG. 1: Relative permittivity e(i£) of various materials as 
a function of imaginary frequency £ (in units of c//xm) 
or "Matsubara" temperature T — h£ > /2nk-B- Doped sili- 
con corresponds (bottom to top) to doping density pd — 
{1,3,5,10,10 2 } x 10 16 , modeled via an empirical Drude 
model 15 , as is gold [16]. Water, polystyrene, ethanol, teflon, 
and lithium niobate are all modeled via standard Lorentz- 
oscillator models 1171. 



as a function of imaginary frequency z£ (bottom axis) or 
"Matsubara temperature" T = h£/27rk (top axis). Fig- 
ure|2]shows the equilibrium separation d c (in units of /im) 
as a function of temperature T G (0, 400) K (in Kelvin) 
for some of the material combinations (solid/dashed lines 
correspond to stable/unstable equilibria), and demon- 
strates various degrees of T sensitivity. The previously 
studied [9] material combination of teflon/ethanol/silicon 
(data not shown) shows very little T-dependence: d c 
varies < 1% over 400 K. More dramatic behavior is 
obtained for lithium niobate (LiNb03) or doped silicon 
(doping density p d = {1, 10, 100, 500, 10 3 } x 10 17 ), whose 
low-£ e crossings with ethanol lead to the desired os- 
cillatory /(£) in (fil: the stable-equilibrium separation 

(s) 

d y c '' for both cases decreases by > 2 pm over 400 K, and 
^d ( c s) w -2 nm/K near T = 300 K (where d { c s) - 300- 

700 nm). The sign of ^d^ comes from the increas- 
ing domination of the repulsive small-£ (large-separation) 
contributions to /(£) as T increases. Varying the sili- 
con doping density dramatically changes the T depen- 
dence because it tunes the low-£ silicon/ethanol e cross- 
ing in Fig. [T] Doped-silicon exhibits another interest- 
ing behavior: d c disappears at a critical temperature 
T c (determined by pd) due to a bifurcation between the 
stable/unstable equilibria. (T c can be tuned not only 
by changing pd but also by changing the LiNb03 layer- 
thickness, here ~ 50nm.) Experimentally, such a bifur- 
cation yields an irreversible transition from suspension 
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FIG. 2: Equilibrium separation d c (in units of /xm) as a func- 
tion of temperature T (in Kelvin), for a geometry consisting 
of fluid-separated semi-infinite slabs (no gravity). The var- 
ious curves correspond to d c for various material combina- 
tions. Solid/dashed lines correspond to stable/unstable equi- 
libria, and shaded regions are T where ethanol is non-liquid 
at 1 atm [18]. Doped-silicon is plotted for various doping 
densities p d = {1, 10, 100, 500, 10 3 } x 10 17 . 

(T < T c ) to stiction (T > T c ). Figure [2] also shows 
a small sample of the many other material possibilities. 
The shaded regions in Fig. [2] correspond to T above the 
boiling point (320 K) or below the freezing point (159 K) 
of ethanol at 1 atm [T8] . 

The inclusion of gravity /buoyancy introduces another 
force into the system and leads to the possibility of 
additional phenomena, such as additional stable equi- 
libria due to gravity/ Casimir competition [9 j. For ex- 
ample, Fig. [3] shows the equilibrium separations d c 
of a polystyrene (PS) slab of thickness h in ethanol 
above a semi-infinite doped-silicon slab (pd = 1.1 x 
10 15 ), including gravity (mass density p PS - Methanol = 
0.264 g/cm 3 [18 ). As in Fig. § d c varies dramatically 
with T: ^d c « 1.2 nm/K near T = 300 K. Gravity be- 
comes increasingly important as h grows: compared to 
h = (leftmost line), it creates an additional stable equi- 
librium (solid lines) at large d c (hundreds of nm) where 
the downward gravity dominates. With gravity, there are 
three stable/unstable bifurcations instead of two, leading 
to three critical temperatures where qualitative transi- 
tions occur: T g refers to the temperature of the topmost 
bifurcation, created by gravity, and the other two tem- 
peratures are labeled T x (« 100 K) and T 2 (« 180 K). If 
T g <T\ (h < 40 nm), there exists an irreversible transi- 
tion from suspension to stiction as T is decreased below 
T g . If Ti < Tg < T2 , there are two irreversible transitions 
from suspension to suspension (smaller d c ) to stiction as 
T is lowered from T > T g to T < T\ starting in the large- 




temperature T (K) 

FIG. 3: Equilibrium position d c (in units of /im) of a 
semi-infinite polystyrene (PS) slab immersed in ethanol 
(shaded T — non-liquid) and suspended against gravity by 
a repulsive Casimir force exerted by a doped-silicon (Si) slab. 
The solid/dashed lines correspond to stable/unstable d c , and 
each color represents a different value of PS slab-thickness 
h (in units of /im). The inset shows the magnitude of the 
total energy Uri^d) (in units of k^T) as a function of d for 
h = 150 nm, at various T. 

d c equilibrium. Finally, when T g — >• T2 (h ~ 300 nm) the 
two stable equilibria merge and only the T\ bifurcation 
remains. Perhaps most interestingly, when this merge oc- 
curs the slope -§^d c can be made arbitrarily large but fi- 
nite, corresponding to an arbitrarily large (but reversible) 
temperature dependence. For example, Ad c « 130 nm 
for a small change AT « 5 K around T2, for h = 300 nm. 

In a real experiment, the situation is further com- 
plicated by Brownian motion, which will cause the 
separation to fluctuate around stable equilibria and 
will also lead to random transitions between equilib- 
ria [20]. In the example of Fig. [3J the attractive in- 
teraction at small separations means that there is a 
nonzero probability that the slabs will fluctuate past the 
unstable-equilibrium energy barrier A/7 T into stiction, 
but the rate of such a transition decreases proportional 
to exp(— A/7 t //cbT) [20] — here, assuming a 50 x 50 /mi 2 
PS slab, ALT T //c B T « 10 4 , so the stiction rate is negli- 
gible. The energy landscape Ut(cI) /k&T is plotted for 
several cases in the inset to Fig. [3j the general prediction 
of experimental observations involves a viscosity-damped 
Langevin process [20] that is beyond the scope of this 
paper to model, but by choosing T one can make the 
potential barrier between the two stable equilibria arbi- 
trarily small and therefore should be able to reach an 
experimental regime in which "hopping" is observable. 

Alternatively, we consider a simpler example system 
with only a single stable equilibrium and a single degree 
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FIG. 4: Average separation (d) (circles) and equilibrum sep- 
aration d c (red line), in units of /im, vs. temperature T (in 
Kelvin) , for a geometry consisting of a fluid-separated hollow 
PS sphere of inner/outer radius r/R = 3.2/5/im suspended 
in ethanol against gravity by a doped-silicon slab and subject 
to Brownian motion. Shaded region indicates where sphere is 
found with 95% probability. The thin black line is the average 
(d)i4o if the Casimir energy at 140 K is used instead of the 
true temperature-dependent energy landscape. 



Experimentally, measuring hundreds of nm changes 
in separation over tens or hundreds of Kelvins appears 
very feasible, perhaps even easier than traditional mea- 
surements of Casimir forces. (In a fluid, static-charge 
effects can be neutralized by dissolving electrolytes in 
the fluid [14], which also have the added benefit of sig- 
nificantly reducting/increasing the freezing/boiling point 
of the fluid [T5].) Such temperature-dependent suspen- 
sions may even have practical applications in microflu- 
idics. We believe that the examples shown in this let- 
ter only scratch the surface of the possible tempera- 
ture/dispersion effects that can be obtained in Casimir- 
suspension systems. Not only are there many other pos- 
sible materials and geometries to explore in the fluid 
context (along with more detailed calculation of the 
Brownian dynamics), and by no means are the effects 
shown here the maximum possible, but similar princi- 
ples should apply in other systems exhibiting competing 
attractive/repulsive Casimir- force contributions. 
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D-0004, by US DOE Grant No. DE-FG02-97ER25308, 
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of freedom: a hollow PS sphere (experimentally avail- 
able at similar scales [21]), filled with ethanol, of in- 
ner/outer radius r/R = 3.2/5/im suspended in ethanol 
above a doped-silicon (pd = 1.1 x 10 15 ) substrate, shown 
on the inset of Fig. [4j (To compute the Casimir energy 
in this system, we employ a simple PFA approximation 
that is sufficiently accurate for our purpose. Here, for 
d ~ 500 nm, the exact energy is « 85% of the PFA 
energy.) For this example, in Fig. [4] we plot the mean 
surface-surface separation (d) ~ J dz z exp[/7 T (2;) / k^T] 
(determined only by the energy landscape and the Boltz- 
mann distribution [20]), corresponding to an experiment 
averaging d over a long time, along with a confidence in- 
terval (shaded region) indicating the range of d where the 
particle is found with 95% probability. The sphere ex- 
periences an attractive interaction at small separations, 
but again we find that the unstable-equilibrium energy 
barrier is sufficiently large (MJ/k&T « 50) to prevent 
stiction for T near 300 K. As T varies, two factors af- 
fect (d): the T-dependence of the Casimir energy Ut(z), 
and the explicit k^T in the Boltzmann factor. To dis- 
tinguish these two effects, we also plot (thin black line) 
(d)i4o ~ / dzzexp[U U o(z)/k B T] where the T = 140 K 
(bifurcation point) Casimir energy is used at all tempera- 
tures. Comparing (d) with (d)i4o, it is evident that most 
of the positive-slope T dependence of (d) (« 0.8 nm/K 
around 300 K) is due to /7 T , and therefore (d) offers a 
direct measure of the Casimir-energy T dependence. 
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